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ABSTRACT 

Since many or most galaxies have central massive black holes (BHs), mergers of galax- 
ies can form massive binary black holes (BBHs). In this paper, we study the evolution 
of massive BBHs in realistic galaxy models, using a generalization of techniques used 
to study tidal disruption rates around massive BHs. The evolution of BBHs depends 
on BH mass ratio and host galaxy type. BBHs with very low mass ratios (say, ^ 0.001) 
are hardly ever formed by mergers of galaxies because the dynamical friction timescale 
is too long for the smaller BH to sink into the galactic center within a Hubble time. 
BBHs with moderate mass ratios are most likely to form and survive in spherical or 
nearly spherical galaxies and in high-luminosity or high-dispersion galaxies; they are 
most likely to have merged in low-dispersion galaxies (line-of-sight velocity dispersion 
^ 90kms~^) or in highly flattened or triaxial galaxies. 

The semimajor axes and orbital periods of surviving BBHs are generally in the 
range 10~^-10pc and 10-10^ yr; and they are larger in high-dispersion galaxies than 
in low-dispersion galaxies, larger in nearly spherical galaxies than in highly flattened 
or triaxial galaxies, and larger for BBHs with equal masses than for BBHs with un- 
equal masses. The orbital velocities of surviving BBHs are generally in the range 
10^-10'*kms~^. The methods of detecting surviving BBHs are also discussed. 

If no evidence of BBHs is found in AGNs, this may be either because gas plays a 
major role in BBH orbital decay or because nuclear activity switches on soon after a 
galaxy merger, and ends before the smaller BH has had time to spiral to the center of 
the galaxy. 

Key words: black hole physics - galaxies: evolution - galaxies: interactions - galaxies: 
kinematics and dynamics - galaxies: nuclei 



1 INTRODUCTION 



It is believed that many or most galaxies house massive black holes (BHs) at their centers (e.g. Magorrian et al. 1998). 
Mergers of galaxies, which are a central part of the galaxy formation process in the hierarchical structure formation picture, 
will inevitably form binary black holes (BBHs) if every galaxy has a central BH and the BH inspiral time is less than a Hubble 
time. The questions of whether the binaries merge and how long they survive before the merger are relevant to resolving a 
number of problems in extragalactic astronomy, such as the detection of gravitational waves (e.g. Folkner 1998) , the bending 
and precession of radio jets from active galactic nuclei (AGNs) (Begelman, Blandford & Rees 1980), the shapes of Fe Ka 
emission line profiles from AGNs (Yu & Lu 2001) etc. On the other hand, so far there is not much observational evidence that 
BBHs exist, which might be expected if (i) the lifetime of the BBH is much shorter than the Hubble time, (ii) the hierarchical 
model of galaxy formation is incorrect, or (iii) the BHs are efliciently ejected from the galaxy by three-body interactions, or 
other processes. The purpose of this paper is to study the lifetime of BBHs in realistic galaxy models. 

Consider two massive BHs with masses mi > 7712 in the core of a galaxy. The merger of the BBH mainly involves four 



processes (Begelman, Blandford & Rees 1980). First, the dynamical friction stage: each BH sinks independently towards the 
center of the common gravitational potential on the Chandrasekhar dynamical friction timescale: 
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(1) 



where the core is assumed to have a one-dimensional velocity dispersion (Jc, a core radius and to contain A'^ stars. Second, 
the non-hard binary stage^ the two BHs eventually form a bound binary system that continues to lose energy by dynamical 
friction. As the binary shrinks, the effectiveness of dynamical friction slowly declines with the increase of the BH velocity and 
the shortening of its orbital period, and the three-body interactions between the BBH and the stars passing in their vicinity 
gradually become the dominant factor that makes the BBH lose energy. Third is the hard binary stage: when 
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(G; gravitational constant), the BBH becomes hard ( Cjuinlan 1996| ). Hard BBHs lose energy mainly by interacting with stars 
passing in their vicinity, most of which will be expelled from the BBH with an energy gain after one or more encounters with 



it. In such expulsions the average relative change in the BBH energy, AE/E, is independent of E (Heggie 1975). The final 



stage is the gravitational radiation stage: as the BBH continues to harden, its semimajor axis decreases to the point at which 
gravitational radiation becomes the dominant dissipative force. A BBH on a circular orbit will then merge within the time 
( |Peters 196^) : 
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where /ii2 = mi 7712 /(mi + m2) and M, = mi + m2 are the reduced mass and total mass of the binary, respectively, and 
m2 < mi. If m2/mi ^ 1 (hence, Oh in equation |2| is small), gravitational radiation can become the dominant dissipative force 
before the hard binary stage begins (see Fig. |^ or § [3.3.4 below) . 

The largest uncertainty in the BBH lifetime comes during the second (non-hard binary) stage or the third (hard binary) 
stage, either of which can be the slowest of the above four processes. The time spent in the second or third stage may vary by 
several orders of magnitude depending on whether or not the low-angular momentum core stars that interact most strongly 
with the BBH are depleted before the BBH decays to the gravitational radiation stage. 

In this paper, we will estimate BBH evolution timescales in realistic galaxy models obtained from observations of the 
central regions of nearby galaxies. To determine the BBH evolution timescale at the slowest stage, the crucial thing is to get 
the BBH hardening timescale associated with interactions with stars in the hard binary stage (see equation ^ below), which 
is longer than the hardening timescale in the non-hard binary stage. The BBH hardening timescale in the hard binary stage 
gives the slowest evolution timescale or its upper limit no matter whether the slowest stage is the hard binary stage or the 
non-hard binary stage. So, our analysis starts from the beginning of the third of the above stages, when the BBH becomes 
hard. The theory we shall develop in § ^ and § ^ has been applied to study tidal disruption rates around massive BHs, the 
evolut ion of galactic nuclei, etc. and was first developed to describe the e volution of globular clusters containing massive black 
holes (Lightman 



Shapiro 1977 



Cohn & Kulsrud 1978 



Shapiro 1985). A simple analysis of the timescales in the first two 
stages is given in § ^ We will only consider purely stellar dynamical processes, ignoring the uncertain, but usually small, 
contribution from gas. We will also ignore the response of the core structure to the BBH evolution and assume that the stellar 
system is in a fixed steady state. The justification for this assumption is that the total stellar mass removed from the galactic 
core during the BBH hardening is generally less than the core mass (see equations ^ ^ or Figures ^ ^ below) . We present 
the results of the evolution of BBHs in realistic galaxy models in § ^ and the estimated observational properties of surviving 
BBHs in § ^. Finally, discussion and conclusions are given in § ^. 



2 LOSS REGION 

In this section, we introduce the concept of "loss region" for stellar systems with BBHs; more specifically, we shall use the 
terms "loss cone" for spherical systems and "loss wedge" for axisymmetric systems. 



2.1 Spherical systems 

In a stellar system with a single stellar mass m, and radius r«, the distribution function (DF) /(x,v) is defined so that 
/(x, v) d'^xd'^v is the number of stars within a phase-space volume d^xd'^v of (x, v). We will consider the generalization to a 



distribution of masses and radii in § 5.6. By Jeans's theorem, in a spherical system the DF depends on (x, v) only through 



* This stage is not called the "soft" binary stage because as Quinlan (1996) pointed out, the term "soft" is best restricted to the familiar 
sense given by Heggie's law, "soft binaries grow softer" (i.e. gain energy and semimajor axis) and a massive BBH in a galaxy core never 
grows softer. 
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the orbital binding energy per unit mass £ — tpi^) ~ ''■nd angular momentum per unit mass J = |x x v|, where ip{r) is 
related to the gravitational potential $(r) through i/'(r) = — <l>(r). If a BH has mass A/,, stars that pass within a distance 
n ~ {M,/mtY^^rt of the BH will be tidally disrupted if rt ^ rs = 2GM,/c? and swallowed whole if rt ^ r^, where c is the 
speed of light, and is the Schwarzschild radius. The region in the (specific energy, specific angular momentum) phase space 
where stars can be tidally disrupted at pericenter is called the "loss cone" (cf. Frank & Rees 1976), and given by: 

<jU£) = '2r\[4>{n)^£]^2GM.n {£<.GM./rt). (4) 

Now, consider a hard BBH with the masses mi > m2 ^ m* and semimajor axis a ^ ay, (see equation ^). The radius of 
the sphere of influence of the BH, an, is defined implicitly in terms of the intrinsic one-dimensional velocity dispersion of the 
galaxy oir) through 

cr^(aH) = GM./an = o-H- (5) 

Note that oh is different from ah = (m2/4M.)aH which represents the semimajor axis at which the BBH becomes hard. 
Given a total BBH mass M., ah depends on the mass of the smaller BH m2, while an does not. We set the core velocity 
dispersion tJc = o"h since in the core, the one-dimensional velocity dispersion (j{r) generally varies only slowly outside an 
[e.g. the "77" models of the density distribution for 1 < < 3 in Dehnen (1993) and Tremaine et al. (1994)]. At large radii 
(r ^ o), the BBH acts as a central point with mass M, = mi + m2 and the stars move in the potential of the stars and 
the central mass point Af, , with no systematic loss or gain of energy. The stars interacting most strongly with the BBH have 
low angular momenta with pericenters /a a (/a is a dimensionless factor ~ 1). They may either lose or gain energy after 
their first encounters with the hard BBH, but eventually most of the stars will be expelled from the potential of BBH with 
an energy gain (Quinlan 1996). As will be seen from equation (^^ below, the energy gain of a star is usually large enough 



for it to escape from the core when the BBH is hard, especially when a ^ ah. By replacing rt with faa in equation (|J), the 
loss cone in the BBH system is given as: 

< JU^JcO) ~ 2GM.faa = 2G(mi + m2)faa {£ < GM./fao)- (6) 

thus Jic will decrease as the BBH hardens. The BBH interacts most strongly with the core stars in the loss cone, and the 
stars escaping from the core with energy gain are considered as being removed from the loss cone. Q The BBH energy loss 
rate is determined by the rate of removal of stars from the loss cone (the "clearing rate"). At first, the depletion of the initial 
population of stars in the loss cone determines the clearing rate. As stars are removed from the loss cone, new stars are 



199e 



scattered into the loss cone by two-body relaxation (Binney & Tremaine 1987) and resonant relaxation (Ranch & Tremaine 



Ranch & Ingalls 1998), and eventually the clearing rate reaches a steady state controlled by the balance between the 
loss rate and the rate at which stars refill the loss cone. If the rms angular momentum transferred to or from the stars per 
orbital period is larger than Jic, then the stars will refill the loss cone as fast as it is depleted and the loss cone remains 
full; otherwise, the stars will slowly diffuse into the loss cone and the loss cone remains nearly empty. Thus at large radii, 
the loss cone is full ( "pinhole" regime) and at small radii, the loss cone is nearly empty ( "diffusion" regime) [the terms come 
from Lightman & Shapiro (1977)]. The transition radius or binding energy between the two regimes is denoted by ric or £ic 
[^(ric) = £ic\. If most of the stars contrib uting to the BBH orbital decay have energy £ <^ GM, / faO., which is generally a 



good approximation as will be seen from § 5.3.4 , the rate of refilling the loss cone in a BBH system can be determined from 
existing tidal disruption calculations, simply by replacing the tidal disruption radius rt with faa. The refilling rate caused by 
two-body relaxation can be obtained by solving the steady-state Fokker-Planck equation (e.g. equation 17 in Magorrian & 
Tremaine 1999, hereafter MT). The loss cone in the BBH system can be approximated as static even though the loss cone 
shrinks slowly as the BBH orbit decays because the timescale of the stellar diffusion into the loss cone is generally shorter 



than the BBH hardening timescale (see § 5.3.4 below). Resonant relaxation is effective only within a radius enclosing a mass 



0.1 Af, of stars and does not contribute significantly to the total refilling rate; therefore, it will be neglected (see MT or 



5.3.4 below). 



2.2 Axisymmetric systems and triaxial systems 

In an axisymmetric galaxy, the DP involves three integrals: /(£,Jz, J3) where Jz is the component of angular momentum 
along the symmetry (z-)axis and J3 is a third integral. If the axisymmetric galaxy is nearly spherical, the third integral J3 may 
be approximated by the angular momentum J (Binney & Tremaine 1987). In axisymmetric galaxies, there exist centrophilic 
orbits such as box orbits, which pass arbitrarily close to the center and have low angular momentum, as well as centrophobic 



t The assumption that stars escaping from the core are removed from the loss cone is also made by Begelman et al. (1980), Quinlan 
(1996), and others. This assumption is plausible since such stars have much longer periods and are much more susceptible to external 
torques than core stars. However, in some circumstances even a small fraction of stars returning to the loss cone could significantly 
enhance the BBH decay rate. Therefore this assumption deserves further investigation. 
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orbits such as loop orbits, which avoid the center and have high angular momentum. Here, we introduce Js to mark the 
transition from centrophilic (J Js) to centrophobic (J <: Js) orbits. Stars on centrophilic orbits with \ Jz\ < Jic can precess 
into the loss cone, while those on centrophobic orbits with J > Jic cannot. Hence, in an axisymmetric galaxy, the loss cone is 
replaced by a loss wedge with \Jz\ < Jic (MT). The BBH energy loss rate is determined initially by the clearing rate from the 
loss wedge and in the steady state by the rates at which two-body relaxation and tidal forces refill the loss wedge. 

Likewise, in a triaxial galaxy, we may expect that there also exists some characteristic angular momentum Js which marks 
the transition from centrophilic to centrophobic orbits, and most of the stars with J < Js can precess into the loss cone. We 
shall call the region J < Js the "loss region" in triaxial galaxies. Thus, the BBH energy loss rate is determined initially by 
the clearing rate from the loss region and in the steady state by the rates at which two-body relaxation and tidal forces refill 
the loss region. 



3 EVOLUTION TIMESCALES AFTER BBHS BECOME HARD 



The stars contributing to the orbital decay of a hard BBH come from the loss region. Assuming that F{£,a)d£ is the clearing 
rate from the loss region for stars with energy £ ^ £ + d£, the energy loss rate of the BBH is given by: 

dE , 



dt 



(a) 



A£F{£,a)d£, 



(7) 



where A£ is the average specific energy change of the stars escaping from the BBH during an interaction with the BBH, and 
the semimajor axis of the BBH is a. The quantity A£ is only weakly dependent on £ for hard binaries. Hence we assume 



A£ 



~KGfii2/a, 



(8) 



where fii2 is the reduced mass of the BBH and K is a constant which will be determined from Quinlan (1996) in § 3.1. Thus, 
after the BBH becomes hard, its hardening time due to stellar interactions is given as: 

M. 1 



th{a) = 



a 




E 


a 




E 



2Km-^ J F{£,a)d£' 



The timescale associated with gravitational radiation is given by (Peters 1964) 



gr(a) = 



64 GVi2M2(l + 73/24e2 + 37/966-1) ' 



(9) 



(10) 



where e is the BBH orbital eccentricity. The merger timescale in equation ^ is related to the timescale in equation (^o|) by 
tmcrgc(a) = Jg tgi{a, e = 0) da/ a. 

As will be seen in § 3.1, the hardening timescale th{a) of a hard BBH is independent of its orbital eccentricity e. The 
gravitational radiation timescale in equation ( p^ depends weakly on e when e is small [e.g. tgr(a, c = 0.3)/fgr(a, e = 0) ~ 0.6]; 
but for large e, it significantly decreases [e.g. tgi{a,e — 0.8)/fgr(a,e = 0) ~ 0.01]. In this paper, we assume that the initial 
eccentricity when the BBH becomes hard is small (say, e ^ 0.3). In this case, in the hard binary stage, the BBH eccentricity 
hardly grows as the BBH hardens (Quinlan 1996); in the gravitational radiation stage, the eccentricity decays exponentially. 
So in our calculations we will always set e = 0. 

Combining equations 1^ and (^o|), we find that the evolution timescale of the hard BBH tevoi(a) is given by: 

1 _ 1 ^ 1 ^^^^ 
tevoi(a) th{a) tgr{a) 

and gravitational radiation becomes dominant when a < Ogr, where agr is defined by fh(agr) = igr(agr). This marks the 
transition between the hard binary and gravitational radiation stage. The total mass of the stars interacting with the BBH 
during the hard binary stage (ogr < a < ah), which can be derived from equation (W, is given by: 



Ml" 



F{£,a)d£dt=^H^: 



(12) 



3.1 K 

Quinlan (1996) studied the dynamical evolution of massive BBHs by scattering experiments using the restricted three-body 
approximation. Consider a hard BBH in a galaxy core with uniform density p and Maxwellian velocity dispersion (Tc. The 
stars in the core start at a radius r ^ a from the BBH, with velocity v. Assume that Vi and Vt are the radial and tangential 
components of the vector v. At large radii, the low angular-momentum stars that contribute to the orbital decay of the hard 
BBH are in nearly radial motion and the number flux into the loss cone is: 
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(13) 



1 



■ exp(--!;^/2cr^). 



m. (27ra|)3/2 

With Ji^lfao) ~ 2GM, faa at large radii, we have from equation (^) the hardening timescale of the BBH: 
, , . a 1 (Jc 

a 4:V2^Kfa Gap 
In Quinlan's simulation, the hardening time is given as 

(Tc \ /10^Mq/pc3\ /0.1pc\ /16^ 



th = 



'^'^ = 2.8 X 10'' yr 



GpaH 



200kms-i 



(14) 



(15) 



(16) 



where for a hard BBH, _ff is a constant (~ 16) independent of the BH mass ratio and orbital eccentricity. Combining equations 
(^ and (|l|), we have the constant Kfa = 0.097// = 1.56. In our calculations, we will always set K = 1.56 and fa = 1. Thus, 
combining equations ^ and (^, the average specific energy change is given by 

2mi \ f o,h\ ( 3(7^ " 



a \ 1.56/ Vmi + m2 

which is generally large enough for the star to escape from the core when the BBH becomes hard, especially when a ^ ah 



(17) 



3.2 F{£,a) 

Most of this subsection is a summary and generalization of the tidal disruption rates calculated in MT. 



3.2.1 Spherical galaxies 

In a spherical galaxy with a DF f{£,.J^), when the BBH first becomes hard (the time t is set to 0), the total stellar mass in 
the loss cone is given by: 

Mic(ah) ~ m. y A-K^T){£)f{£,f = Q)JI{£ Jaai,)P{£) A£ < ri{£) < 1, (18) 

where P{£) is the radial period of an orbit with energy £ and zero angular momentum and ri[£) is a dimensionless factor. If 
'rj[£) = 1 (or 0) for all £, the loss cone is full (or empty). The clearing rate per unit energy at time t is first mainly controlled 
by the draining rate from the loss cone: 

^drain.p ^. ^ [ ^^^viniiS, = 0) jf,(£:, /.a) A£ for t < P{£) 

^ ' ' \ for f > P{£) ' ^ ' 

After time P{£), when the loss cone of stars with energy £ is depleted, the loss cone is refilled by two-body relaxation. The 
steady-state diffusion rate of stars into the loss cone is given by the steady-state solution of the Fokker-Planck equation (13) 
in MT, which is a generalization of the equation in Cohn & Kulsrud (1978) to a non-Keplerian potential: 

F-i£,a)d£^ /^Zf^f (20) 

where 

/^^(f) = A7v^fi£)J^{£)p{£)Pin, (21) 

R (F f r,)= -^iciSJaa) ^ / eM-liSJaa)] for q{£Jaa) > 1 

Uo(C,jaa)- ^2^^^ \exp[-0.186g(£:,/„a)- 0.824y^(?JU)l for ^(f, /„a) < 1 ' ^^^> 

q{£,faa) = m£)Pi£)j!i£)/jU£,faa), (23) 

'-'•^'-f'^M '' = ^- 

Jc(f) is the specific angular momentum of a circular orbit at energy £, f{£) is the "isotropized" DF defined by f{£) = 

j2 

Jq ° f{£) J ) d>/ / Jc (f ), and r- are the apocenter and pericenter of loss-cone orbits, (A«t ) is the diffusion coefficient for 
tangential velocity (see equation 8-64 of Binney & Tremaine 1987), and p[£) is the orbit-averaged diffusion coefficient. Note 
that the stellar encounters responsible for the velocity diffusion can be treated as elastic encounters only outside the "collision 
radius" : 
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rcoii ==; 7 X lO"'(A/./m.)(r./R0)cm (25) 



at which the velocity dispersion ~ ^JUm^Jt is comparable with the escape velocity from typical stars (Frank & Rees 1976). 
When r ^ r-coii, two stars cannot deflect each other's velocities through a large angle without coming so close that they 
actually collide. Generally, most stars have Ro{£,faa) ^ 1, and the isotropized DF is close to the real stellar DF and is 
related to the stellar density p{r) by 



p(r)«4™. / f{£)^2[<5>ir)^£]d£. (26) 



The transition between pinhole {£ < £ic) and diffusion {£ > £ic) regimes is at q{£ic,faa) — — ln[Ji^(fic, /aa)/Jc (^■)1- From 
equations (po|)-(p^ it can be seen that in the diffusion regime the flux into the loss cone F^'^{£,a) is insensitive to a since it 
depends only logarithmically on a. 

Considering both the draining and the refilling of the loss cone, at a certain time t, the clearing rate is given as: 

F(f,a;t) =max[F'''^"'"(£:,a;t),F'"(£:,a)]. (27) 



3.2.2 Axisymmetric galaxies 

In spherical systems, stars are all on centrophobic loop orbits. In axisymmetric galaxies, there also exist centrophilic orbits 
which can enhance the clearing rates. The dynamics of eccentric orbits in non-spherical potentials can be studied by a simple 
symplectic map constructed by Touma & Tremaine (1997). The mapping models the evolution of eccentric orbits as a two-step 
process: (i) precession of the orientation of the orbit in a non-spherical potential, and (ii) a kick to the angular momentum of 
the orbits at apocenter where the star spends most of its time, and the torques are likely to be strongest. Following MT, we 
use this mapping to study the dynamics of centrophilic orbits with = in axisymmetric galaxies: 

K = Yn- ^esin26l„ 



Jn + l 



= e„ + 5(K) 



K+i = K-iesin2e„+i (28) 

where Y = J±/ ,Jc, J± is the scalar angular momentum {\J±\ — J) which can be positive or negative, 9 is the colatitude at 
each apocenter passage, eJc is the time-integral of the torque over one radial period and g{Y) is the advance in 9 per radial 
period. The parameter e is related to the axis ratio of the potential b and the mass ellipticity e', for example, 

e = \/2to(1 - 6) = V2^e/3 (b < 1) (29) 



in the logarithmic potential (Touma & Tremaine 1997) 



ilog(x?+xi/62). (30) 

Of course, the potential of realistic galaxies with central BHs is not scale-free, and so the precession rate and the mapping in 
equation (^^ are dependent on the energy £. As shown in Figure 4 of MT, in realistic galaxies, the centrophilic orbits are 
regular at large £ and stochastic below some critical energy fcrit. Above fcrit, the peak angular momenta of the regular orbits 
passing the points (F = 0, ^ = 0,7r) in the Y-9 surfaces of section are denoted by Jm(£) and the area of {Y,6) phase space 
covered by the regular orbits with peak angular momenta less than Jm is 4Jn-i{£) / Jc{£) ■ Below fcrit, the areas of {Y, 6) phase 
space occupied by the stochastic orbits are denoted by 2ttJi{£) / Jc{£). Those stars on centrophilic orbits with \ J±\ < Jm (for 
regular orbits) or | J±| < J; (for stochastic orbits) can precess into the loss cone so long as | Jz| < Jic. Here, we use one symbol 
Js{£) to represent 2Jni(£)/vr for regular orbits or Ji{£) for stochastic orbits. 

If the DF of the axisymmetric galaxy is a function of two integrals f{£, Jz), when the BBH first becomes hard, the total 
mass in the loss wedge is given by: 

Mi„(ah)~m.y 'iTv^rj{£)f{£,J, =0)m&x[jU£,faak),2Ji,(£,faai,)Js{£)]Pi£)d£ < r]{£) < I, (31) 

where as in equation (jl^), the dimensionless factor ri{£) = 1 (or 0) for a full (or empty) loss wedge. If Jic{£, faO) is larger than 
2Js{£), we use the clearing rates obtained in the spherical case. After Jic{£, faa) becomes less than 2Ja{£) at a time given by 
t = Ts{£) (t = at a — tth), we will consider the effects of flattening as follows. First, the loss wedge in a two-integral model 
drains at a rate: 



^d.ain,^ tWf ~ ^ 47r2,7(f )/(f , J, = Q)Jl{Ua) A£ for t < T<i™(f ) 

^ ' ' ^ ^0 for t > r'*'^-"(£:) ' 
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where T'^''^^'^{£) satisfies 

j2,f ,0 {£,a;t)d£+ dt {a ) ' ' ' ' F {£,a;t)A£ 

•J\c\Ja."- ) Jmin[P{£),Ta{£)l Jlc Ua" ^ Js j 

~ 87r2r,(£:)/(£:, J, = 0) Jie(/aa) Js(f ) df. (33) 

In the above equation, the right-hand-side term represents the total number of stars (with energy f — > f -|- d£) which can 
precess into the loss cone J < JidfaCi) by tidal forces. If Ts{£) < P{£), the left-hand-side terms in equation ( |33| ) approximately 
represent the numbers of stars (with energy £ £ + A£) removed from the loss wedge before and after J\c{£, fai) — 2Js{£); 
and if r,(£:) > P{£), we have ydrain^^) ^ p^^^y 

After the loss wedge is depleted, the rate at which the loss wedge is refilled is given by a Fokker-Planck analysis in MT: 

F^-i£,a)A£^4.^jU£,f.a)^^ij^)^^, (34) 

where = q/8 (see q in equation and f{£) is the "isotropized" DF defined by f{£) = J^" f{£, Jz) dj^/ Jc{£)- If the galaxy 
is close to spherical, the isotropized DF f{£) is close to the real stellar DF and f{£) is also related to the stellar density by 
equation (^6|). The transition between the "diffusion" and "pinhole" regions occurs at q ~ 2Jc/Js. From equations (^) and 
( p^ it can be seen that in the diffusion regime the flux into the loss wedge F^^{£,a) is independent of a. 

Considering both the draining and refilling of the loss wedge or the loss cone, at a certain time t, the clearing rate is 
given as: 

F(£:,a;t) =max[F'''^"'"(£:,a;t),F'"(£:,a),F''^(£:,a)]. (35) 



3.2.3 Triaxial galaxies 

Similarly, in triaxial galaxies, there will be some characteristic angular momentum ,Jb[£) inside which most orbits are cen- 
trophilic and the BBH decay rate is determined by the clearing rates from the loss regions J{£) < Js{£). Here, we approximate 
the Js{£) obtained in the axisymmetric cases as the characteristic angular momentum Js{£) in triaxial galaxies. Merritt and 
Quinlan (1998) argue that the tria:xiality of galaxies with central BHs decays due to the evolution of stochastic orbits, so that 
most galaxies become axisymmetric at any given radius within a time t = ^ IQ^ local orbital periods. When the BBH 

first becomes hard, we set the time t = and the age of the galactic triaxiality is defined as T^gJ, at this time. We give the 
mass in the loss region at a = Oh as: 

Mir{ai,)~m,J 4n^r]{£)f{£)ina.x[jU£,faai,),j'^i£)]P{£)d£ <??(£■)< 1. (36) 
When < t < Ttii{£) = max[r*'^''"° — Tagc, 0], we consider the clearing rates in triaxial galaxies as follows. If Js{£) < Jic{£), 



the clearing rates follow the same formula as in the spherical case (see § 3.2.1). When Js{£) becomes larger than Jic{£) at a 



time given hy t = Ts(£), the clearing rates from the loss region are first controlled by the draining rates from the loss region: 
^d.ain,^ . ^ Un^7^i£)f(£)Jl{faa)d£ for t < min[T''--(£:), TtH(£-)], 

^ ' ' ^ \0 for f > min[r'''^=''"(£:),rtH(£:)] ^ ' 

where T'*''''" satisfies 

mm[P(£),T,(£)] j'^ ( f „T<1™° 

dt'{a') iiy, ' F'"''''''{£,a';t')d£+ / dt' {a') F"'"'''' {£ ,a' -t') d£ 



Jlcifo.°-') ^miii[F(£),Ta{£)l 



^.AiT\{£)f{£)j''A£)P(£)d£. (38) 



Equation ( |38[ ) has similar meanings to equation (|33|). 

After the loss region is depleted, the refilling of the loss region can be described by equation (^^ or (p^. The clearing 
rate is given as equation (BST). 



When t > rtri(£'), we will assume axisymmetry and follow the calculation in § 3.2.2 



4 EVOLUTION TIMESCALES BEFORE BBHS BECOME HARD 

In this section, we carry out a simple analysis of the evolution timescales before the BBH becomes hard. Assume that a BH of 
mass mi is located at the center of a spherical galaxy with stellar mass density p(r) and one-dimensional velocity dispersion 
(J if). Suppose that a BH m2 (m2 < mi) that used to be at the center of a galaxy is orbiting with a velocity V2 and spiraling 
into the center of the parent galaxy by dynamical friction. Because the inspiraling BH m2 is accompanied by the stars bound 
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to it with total mass M™^ , dynamical friction brings the two BHs together much more rapidly than if m2 is "naked" (e.g. 
Milosavljevic & Merritt 2001). The dynamical friction force on m2 and its accompanying stars Af™^ is given by (Binney & 
Trem4.ine 1987|): 



(m2 + Mr^) 



/dV2\ 

V At ), 



-47r lnAG^(m2 + Mr^) p(r) 



erf(X) 



(39) 



where X — V2/\/2a{r), erf is the error function, and m2 and its bound stars M™^ are assumed to act as a single body. The 
logarithm of the ratio of the maximum and the minimum impact parameter In A is set to unity (which overestimates the 
dynamical friction timescale, but not by more than a small logarithmic factor). Assuming that the orbit of m2 is circular, the 
force (equation is tangential and causes ?7i2 to lose angular momentum per unit mass L at a rate: 



At V df y df ' 

where 

L(r) = rv2 = v^G[mi + M,(r-)]r 

and M-t(r) is the stellar mass in the parent galaxy within a radius r. The dynamical friction timescale is given by: 



r 




dL/dr 


r 


— J. 


di/df 



(40) 



(41) 



(42) 



which can be calculated from equations (^)-(^l|). 

In the potential of the parent galaxy with BH mi, the stars around the BH m2 are tidally truncated at the radius 

1/3 

n 



4^2 - 
\ 4^2 - K- 



1/3 



»m2, 



(43) 



where Q. and k are the circular frequency and epicycle frequency in the parent galaxy. Note that the tidal forces approach 
zero in a nearly homogeneous core, in which k = 2f2. We use this formula to obtain a crude estimate of the ratio of the 
bound stellar mass M™^ to the BH mass m2. If the distribution of the stars around each of the BHs is a singular isothermal 
sphere and their one-dimensional velocity dispersions are cri and (T2, we have fi^ = 2a\/r'^, = 4cri/r2, Af™^ — 2(7f rt/G 
and rt = (a^j \plrj\)r . If the velocity dispersions a\ and a^. follow the same relation with central BH masses mi and m2 as 
equation (pd) (see § 5.1 below), the ratio of the bound stellar mass Af™^ to the BH mass m2 is given by: 



A^r 



m2 



Gm2Cri 



= 300 



2kpc 



/ m2 
Vmi 



10*^ Me 
mi 



Ml 



» m2. 



(44) 



The stellar mass bound to m2 is effective in speeding the inspiraling of m2 only when m2 is at large radii (e.g. ^ 10-100 pc 
in Fig. ^ below). 

When m2 sinks to the radius at which Af, (r) ~ M, (e.g. r ~ l-lOan in Fig. |^ below) , the two BHs form a "bound" 
binary system (this marks the transition from the dynamical friction stage to the non-hard binary stage). (Here, the definition 
of "bound" is somewhat arbitrary; we only intend to mean that the gravitational force on m^. is dominated by mi rather than 
by the stars.) Though the BBH orbital decay in the non-hard binary stage comes from both dynamical friction from distant 
stars and three-body interactions between the BBH and the stars passing the BBH vicinity, we still use equation (^) to 
estimate the BBH evolution timescales in both the dynamical friction stage and the non-hard binary stage because scattering 
experiments with the restricted three-body approximation basically give a similar hardening timescale to equation (^^ [see 
equations 2, 16 and 18 in Quinlan (1996)]. 

Before the BBH becomes bound [Af, (r) <: Af,], the removed mass from the galactic core can be ignored since both the 
mass and binding energy of the smaller BH are less than the total mass and kinetic energy of the core stars interior to the 
BBH orbit; thus the inspiral has negligible efi'ect on the dynamics of the core stars. But, during the non-hard binary stage, 
depletion of the stars that interact with the BBH can be significant in some galaxies. The mass removing rate from the galactic 
core caused by interactions between stars and one or both of the BHs, Jrm(a), is defined by 

1 dAfrm , ... 

^™(") = M:dwTM' 

where Afmi is the stellar mass removed from the galactic core. In § ^, we introduced the concept of loss region (loss cone for 
spherical galaxies and loss wedge for axisymmetric galaxies) for a hard BBH, but equations (^|), (|^) and (|3^) (i.e., the total 
stellar mass in the loss region) can be generalized to the non-hard binary stage to reflect the mass of the stars which may 
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have close encounters with the BBH. Here, the stars which may have close encounters with the BBH can be either the original 
low-angular momentum stellar population or those precessed from high-angular momentum orbits by tidal forces. Before the 
BBH becomes hard, if the mass of the initial stellar population in the loss region is not enough to provide the removed 
stellar mass from the galactic core Mrm (see equation |4^) , the BBH would lose energy mainly by interactions with distant 
stars, which becomes more and more inefficient as the BBH hardens, and by three-body interactions with the low-angular 
momentum stars which are diffused from high-angular momentum orbits by two-body relaxation, which will dominate the 
BBH hardening timescales in the hard binary stage (see equations ^ and ^). Thus, the BBH hardening timescale should 
be higher than those estimated from equations (p9|)-(|42|) in the non-hard binary stage, and smoothly increase to connect the 
hardening timescale in the hard binary stage. Here, we use the linear approximation [for the ln(fh)-ln(a) relation] to obtain 
the BBH hardening timescale at ah < a < adpi 

ln[th(a)] = ln[th(aii)] + ln[tdf (adp)/th(ah)] ln(a/ah)/ ln(adp/ah), (46) 

where a^p is the semimajor axis where the initial population of stars in the loss region are all removed from the loss region 
(which will be determined in § 5.2), th(a;h) is obtained from § ^ and tdf(a) is obtained from equations (^)-(^). Considering 



the possible BBH energy loss caused by gravitational radiation, the BBH evolution timescale tcvoi(a) is given as equation ( |ll| 
and ih(a) = igr(a) may happen before the BBH becomes hard. 



5 RESULTS 



5.1 Galaxy properties 

In the past several years, images from the Hubble Space Telescope (HST) have revealed many details about the central regions 
of nearby galaxies, with a resolution of 0.1", corresponding to distances of several pc or 10^-10^(M,/10* M©) Schwarzschild 
radii for typical target galaxies (Byun et al. 1996). The inner surface brightness profiles of the galaxies are well fitted with a 
five-parameter fitting function - the Nuker law: 



fi-y 



(47) 



The asymptotic logarithmic slope inside rb is —7, the asymptotic outer slope is — /3, and the parameter a parameterizes the 
sharpness of the break. The break radius rb is the point of maximum curvature in log-log coordinates. The "break surface 
brightness" /b is the surface brightness at rb. Elliptical galaxies and spiral bulges (hot galaxies) can be classified into two 
main types according to their inner surface brightness profiles (7 oc when r — > 0). Hot galaxies with a quite shallow inner 
power-law slope (7 ^ 0.3) are classified as "core" galaxies and those with a steep slope (7 <; 0.5) are labelled as "power-law" 
galaxies. The difference between core and power-law profiles has no direct connection with the presence of a disk, whether 
seen edge-on or face-on, although power-law galaxies have disky isophotes (Faber et al. 1997). Core galaxies are luminous 
(Afv < —20.5) with large central BHs and power-law galaxies are faint {Mv > —22) with small central BHs. Early discussions 
of the demography of the central BHs focused on the correlation of BH mass with galaxy luminosity. It is now believed that 
BH mass is correlated more tightly with galactic velocity dispersion, which suggests a causal connection between the formation 



and evolution of the BH and the bulge (Gebhardt et al. 2000; Ferrarese & Merritt 2000) 



We shall apply the theoretical results of the preceding sections to estimate what the evolution timescales of the BHs in 
the centers of nearby galaxies would be if they were binary. 

We use the HST sample of hot galaxies compiled in the paper by Faber et al. (1997). Among the 61 galaxies in that 
paper, 41 galaxies (Table 1) have values for the Nuker law parameters, the stellar mass-to-light ratio T (constant for each 
galaxy) and the half-light radius, where T is determined by normalizing to the central velocity dispersion based on spherical 
and isotropic models fitted to the Nuker- law profile. The isotropized DF f{£) of the galaxies (see equation Effl can be obtained 



by the Eddington formula (Binney & Tremaine 1987) 
1 d 



where 



^/87^2 d£ 



T 

p(r) = Tj(r-) = 



1 dp d-i/) 
m» Alp y/ijj — E ' 



dJ di? 



The total mass of the central BHs is estimated by: 

/ „ \3.75{±0.3) 

M. = 1.2(±0.2)xlO«Me(^5ol^ 



(48) 



(49) 



(50) 
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where CTc is the luminosity- weighted hne-of-sight velocity dispersion inside the half-light radius (Gebhardt et al. 2000). The BH 
mass-velocity dispersion relation due to Merritt & Ferrarese (2001) uses central velocity dispersion and the velocity dispersion 
has a different exponent, 4.72(±0.36). In this paper, we use equation ( |5o| ) to estimate the central BH mass. The results will 
not be affected much if we equate the central velocity dispersion to CTc and use the different exponent 4.72. Of these 41 galaxies 
in Table 1, 9 DFs obtained from equation ( ji^ ) are negative, perhaps because the spherical and isotropic dynamical model is 
not appropriate or because the Nuker law is a poor fit, and these are deleted from our sample. All of these 9 galaxies are core 
galaxies with shallow central cusps (0^7;$ 0.04). Generally, CTc is determined by the gravitational potential of stars and 
affected little by the gravitational potential of the central BH, but that is not true for two galaxies with 7 > 1, which are also 
deleted from our sample. Thus, a total of 30 galaxies is left in our study. 

5.2 Total stellar mass in the full loss region 

Before studying the details of the BBH evolution timescales, we will first generalize equations (jl^) and ( ^l| ) to the non-hard 
binary stage to obtain the total stellar mass in the full loss cone/wedge (for spherical/axisymmetric galaxies) to see if the 
total stellar mass in the full loss cone/wedge is large enough to provide the removed stellar mass from galactic core before the 
BBH becomes hard (Fig. |). 

Quinlan's simulation (1996) gives the mass removing rate from the galactic core caused by interactions between a BBH 
and stars (Jrm, see equatio n ^5|) as a function of CTc/wc (or a/oh) for five different BH mass ratios (m2/mi=l/256, 1/64, 
1/16, 1/4, 1), where «c = ^GM.ja. We can roughly obtain the removing rates for other BH mass ratios by interpolating the 
rates of those five mass ratios. As shown in Figure |l| (or see Fig. 5 in Quinlan 1996) , the removing rate Jrm is small at large 
semimajor axes (e.g. Jrm 10~^ at a/an — 0.3) and increases steeply as the BBH hardens, and we may approximate equation 
( p5[ ) as Mrm(a) — AJ. Jrm(a). The removing rate is about 0.03-0.1 when the BBH becomes hard, and then increases towards 
a constant in the range 0.1-1 as the BBH hardens. This is roughly consistent with the factor 1/(2J(') ~ 0.3 (Kja = 1.56, 
/a ~ 1) in equation (^ or the energy gain in equation ( |l^ ) (i.e., soon after the BBH becomes hard, especially when a <^ ah, 
most stars having close encounters with the BBH will be removed from the galactic core). 

Figure Q(a) shows the total stellar mass in the full loss cone Mic(a) and the removing rate as a function of a if the galaxies 
in Table 1 are spherical. As seen from Figure ^(a), with the decrease of a, Mic(a) decreases and the mass removing rate from 
the galactic core Jrm(a) increases. For the BBH with equal masses {m2/m\ = 1), we have Mic/Mrm — Afic/Af.Jrm ^ 1 at 
a ^ Oh [i.e., in Fig. |l|a, the solid circle representing the stellar removing rate Jrm(a = flh) for ra^lmx = 1 is located below all 
the Afic(a ^ ah)/M, curves]; so the mass of the removed stars from the galactic core during the non-hard binary stage can 
be ignored and the loss cone is approximately full at a = ah. Compared with the BBH with equal masses, the BBH with low- 
7112 /mi ratio becomes hard at smaller semimajor axis and the removed mass may be larger than the mass in the full loss cone 
before the hard binary stage [i.e., in Fig. solid circles for low-m2/mi ratios are located above some of the AJic(a <: ah)/M, 
curves]. For the BBHs whose loss cone is depleted before the hard binary stage begins (e.g. 20 out of our sample of 30 galaxies 
for ni2/m\ = 1/16), we obtain the depletion semimajor axes adp approximately by setting AJic(adp)/A/, = Jrm(<idp) and we 
will use equation (^6|) to obtain the BBH hardening timescales at Oh < a < ajp in § Note that in Figure ^(a), at a given 
a(Ji O.Ioh), Jrm(ci) decreases with BH mass ratios; and hence ajp also decreases with BH mass ratios. 

Figure 0(b)-(f) shows the total mass in the full loss wedge M\„{a) and the removing rate as a function of a if the galaxies 



in Table 1 are axisymmetric. The total mass in the full loss wedge depends on the degree of flattening e (defined in § 3.2.2), 
which is assumed to be a constant here. The average mass ellipticity of the sample in our study, which can be obtained from 
Table 2 in the paper by Faber et al. (1997), is about 0.26. Their average e is ~ 0.36 if equation ( p^ is applied. As seen from 
Figure ^b)-(f) (e=0.01, 0.05, 0.1, 0.3 and 0.5), Mi^{a) increases with the increase of e. If e is very small (e.g. e < 0.01), the 
loss wedge is approximately full at a = ah only for m2/m\ ~ 1 as in the spherical case. When e ^ 0.1, the loss wedge is 
approximately full at a = ah for 1x12 /m\ ^ 1/256 (i.e., the corresponding solid circle is located below all the M\v,/M, curves 
at a ^ ah in Fig. Qd-f). 

For triaxial galaxies, the maximum mass of the stars which may have close encounters with the BBH is not less than 
the mass in the full loss wedge of axisymmetric galaxies with the same parameter e. In the study below of the evolution of 



BBHs in triaxial galaxies (see § 5.5 below), we will mainly consider BBHs with equal masses. In this case, the removed stellar 
mass during the non-hard binary stage can be ignored and the loss region is approximately full at a = ah as in axisymmetric 
galaxies (cf. Fig. ^-f). For other cases in triaxial galaxies (varying BH mass ratios at a given e), the variation of the BBH 
evolution should follow similar trends to those in axisymmetric galaxies, and we need not consider the maximum mass of the 
stars which may have close encounters with the BBH in triaxial galaxies. 

5.3 Spherical galaxies 

Assuming that the galaxies in Table 1 are spherical, we obtain the BBH evolution timescales as a function of BBH semimajor 
axis (or the separation of the two BHs before they become bound) and the BH mass ratio m2/m\ < 1. The results are shown 
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Figure 1. The total stellar mass in the full loss cone/wedge (for spherical/axisymmetric galaxies) and the removing rate from the 
galactic core as a function of a /an (ajj represents the radius of the sphere of influence of the BH M, , cf. equation The solid 
lines (for core galaxies in Table 1) and the dotted lines (for power-law galaxies in Table 1) represent the total stellar mass in the full 
pes Mic(a) in panel (a) (cf. equation |lq ) and the total stellar mass in the full loss wedges Miw(i) in panels (b)— (f) (cf. equation 
The dot-dashed lines represent the removing rate from the galactic core Jmi (see equation |45| ) for the five different m2 /mi ratios 
{ra2/ra\ = 1/256, 1/64, 1/16, 1/4, 1) which can be obtained from Quinlan's simulation (1996). The vertical dashed lines represent the 
place where the BBH first becomes hard and the corresponding m<2./m\ ratio is labelled at the right side of each dashed line. For each 
given m2/m\ ratio, the crossing of the vertical dashed line (a = a^) and the removing rate curve, Jrm, is labelled as a solid circle. The 
total mass in the full loss region increases with increasing e. The M\^(a\/) and Afl^^,(^^h) also increase with increasing m^/rax ratios. For 
the cases in which the solid or dotted lines are higher than the solid circles at a = ajj, the removed stellar mass fro m galactic cores during 
the non-hard binary stage can be ignored and the loss region is approximately full at a = ajj (see details in § |5.2|). 



in Figure ^, where the BBH evolution curves are composed of two parts: one comes from equations (p|)-(|ll|) and ([l9|)-(|2^) 
in § ^ describing the evolution in the hard binary stage and the gravitational radiation stage (dotted or solid lines), and 
the other from equations (|3^)-(^) or equation ( ^6| ) in § ^ describing the evolution in the dynamical friction stage and the 
non-hard binary stage (long dashed or dot-short dashed lines) . A characteristic feature of all these plots is a pronounced peak 
in the evolution timescale between 0.001 and 10 pc, where dynamical friction and gravitational radiation are both relatively 
ineffective. This is always the slowest evolution stage if we do not consider the possibly longer timescales at very large radii, 
a <: 10* pc. We shall call this stage the "bottleneck". 

To study the effect of different BH mass ratios on the BBH evolution, we will first consider BBHs with equal masses (in 
§ 5.3.1), and then consider those with unequal masses (in § 5.3.2), though actually there is no strict dividing line between 
these cases. 
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Figure 2. BBH evolution timescales for splierical galaxies. The panels show different mass ratios m2/mi and the curves represent 
timescales \a/d\ for the galaxies in Table 1. The long dashed lines (for core galaxies: 7 ^ 0.3) and the dot-short dashed lines (for power- 
law galaxies: 7 ^ 0.5) represent the timescales before the hard binary stage; and the solid lines (for core galaxies) and the dotted lines 
(for power-law galaxies) represent the timescales after the BBHs become hard. In panel (a), the solid circles (for core galaxies) and the 

S;n circles (for power-law galaxies) represent the maximum Brownian motion magnitudes of the BBH centers of mass (rran in equation 
. In panel (b), the solid circles (for core galaxies) and the open circles (for power-law galaxies) mark the boundary between the hard 
binary and gravitational radiation stage or the place where gravitational radiation becomes dominant in the non-hard binary stage (i.e., 
o = ttgr when tji = tgr)- BBHs in power-law galaxies have shorter lifetimes than those in core galaxies. With decreasing mass ratios, the 
evolution timescales at large radii (e.g. a ^ 10^ pc) increase and hence BBHs with very low mass ratios (say, <^ 0.001) cannot be formed 
by mergers of galaxies. The bottleneck stages shift to smaller radii as the mass ratio decreases. In most core galaxies, the bottleneck 
timescales become shorter as the mass ratio decreases; while in most po wer-l aw galaxies, they are affected little by mass ratios. The 



lifetime of BBHs is not affected by their Brownian motion. See details in § 5.3 
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Figure 3. Evolution timescales of BBHs (with equal BH masses) as a function of the BBH separation or the orbit semimajor axes in 
the dynamical friction stage and non-hard binary stage. The solid lines represent core galaxies and the dotted lines represent power-law 
galaxies. At small semi-major axis, the curves are cut off at a = aij. The curves are composed of two groups, which is easy to distinguish 
at large radii (a ~ 10^ — Iff^pc). The lower group {\a/d\ ~ 10*^-10^ yr at a ~ lO" pc) are the evolution timescales including the bound 
stars around m2; and the higher group (|a/d| ~ 5 X 10^"-5 X 10^^ yr at a ~ 10^ pc) are the evolution timescales without including the 
bound stars around m2 (see § Q). At small radii (a — > ah), the two groups have almost the same timescales. At large radii (a ^ 100 pc), 
the evolution timescales decrease with decreasin g a: a t small radii [a <^ 100 pc), the evolution timescales increase for most core galaxies 



and decrease for most power-law galaxies (see § 5.3.1). The evolution timescales at large radii are shorter than the Hubble time 10^" yr 



if and only if we assume that the BHs arc accompanied by the tidally truncated remnant of stars from their original host galaxy. 



5. 3. 1 Evolution of BBHs with equal BH masses 

In this subsection, we will consider the evolution of BBHs with equal BH masses (m2/mi = 1, Fig. ^) for which the removed 
stellar mass from the core duri ng th e non-hard binary stage can be ignored and the loss cone is approximately full at a = Oh 
[ri{£) ~ 1 in equation ^ (see § 5^ or Fig. ^). 



As seen from Figure H(a), the overall evolution of the BBHs has the following trends: starting at large radii (~ lOkpc), 
the evolution timescales first decrease with decreasing radii; then increase at some intermediate radii to reach the bottleneck; 
and finally when the BBH semimajor axis a is small enough 10~^-1 pc), the evolution timescales decrease with decreasing 
a. The bottleneck occurs during the hard binary stage. 

Before the BBH becomes hard, the evolution timescales are always shorter than the Hubble time tnubbic = 10^" yr; these 
relatively short timescales refiect the enhancement of dynamical friction due to the bound stellar mass around m2. Figure ^ 
shows that without including the bound stellar mass around m2, the dynamical friction timescales would be much larger than 
^Hubble at large radii (e.g. ~ 5 x 10^''-5 x 10^^ yr at lOkpc). The dynamical friction timescale decreases with the decrease of 
r, and when r decreases to the place where mi dominates the potential, we have \r/r\ oc v:^^{r)/ p{r) oc r~^'^ /r~'^~^ cx r''""'^ 
(see equations ^^^2|). Thus, at small radii (a ^ 100 pc), with the decrease of r, the evolution timescale continues decreasing 
for power-law galaxies (7 <; 0.5), but increases for core galaxies (7 ^ 0.3). 

As seen from Figure ^(a), after the BBH becomes hard, starting at a = ah, the BBH evolution timescale increases steeply 
with decreasing separation a, then becomes nearly fiat, and then decreases as a* once the BBH reaches the gravitational 
radiation stage. Most of the BBH evolution timescales at the bottleneck are longer than the Hubble time. In these galaxies 
we might expect to find BBHs with semimajor axes in the range 0.01-10 pc. BBHs in the core galaxies generally reach the 
bottleneck at larger semimajor axes (oh ~l-10pc) and need more time to reach the gravitational radiation stage than those 
in the power-law galaxies. At the beginning of the hard binary stage, the BBH orbital decay rate is controlled by the draining 
rate of the stars initially in the loss cone and the BBH evolution timescale increases steeply as the stars in the loss cone are 
depleted. The loss cones are depleted (generally when a ^ O.lan) before the BBHs are close enough for gravitational radiation 
to be effective. Thereafter, the clearing rates are controlled by diffusion of stars into the loss cone. Figure ^ shows that the 
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Figure 4. The semimajor axes a of the BBHs versus the transition radii ric between the pinhole (r > ric) and the diffusion (r < ric) 
regimes in spherical galaxies. The lower ends of the lines are cut off at the Schwarzschild radius rs = 2GM, j r? . The solid and dotted 
lines have the same meaning as those in Figure ^. At a given a, most core galaxies have larger transition radii than power-law galaxies 
because core galaxies generally have more massive BHs and lower central densities. In the hard binary stage, if the loss cone in spherical 
galaxies is depleted, the inner parts of the galaxies are generally in the diffusion regime. 



inner parts of most of the galaxies are in the diffusion regime (r < ric) before the gravitational radiation stage begins (generally 
Ogr/cih ft 0.001 in Figure though the transition radii ric decrease with the shrinking of the loss cone. As mentioned in 
§ |3.2.1 , in the diffusion regime, the flux of stars into the loss cone is insensitive to a, and hence the BBH evolution timescales 
controlled by the stellar diffusion are nearly constant. As seen from Figure ^(a), the bottleneck timescales in core galaxies are 
larger than those in power-law galaxies, which reflects the fact that core galaxies have lower central densities (hence, the stellar 
diffusion rates into the loss cone are slower) and more massive BHs [hence, the average relative change in the BBH energy 
caused by interactions with a star passing the BBH vicinity is smaller {\AE/E\ — 2Kmt/M,)] (see equation^. Gravitational 
radiation becomes dominant at larger semimajor axes in core galaxies, because their BBH bottleneck timescales are longer 
and the gravitational radiation timescale (equation hof) is smaller for large BHs. 
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Figure 5. The ratios of enclosed stellar mass within radius r to BBH total mass. The solid and dotted lines have the same meaning as 
those in Figure ^. The left vertical dashed line represents the semimajor axis aij where the BBH with equal masses first becomes hard. 
The open squares re pre sent the break radii of the surface brightness profile r^, (cf . equation |47| ) ; the open circles represent the collision 
radii rcoii (equation bslli the solid circles represent the boundary between the hard binary and gravitational radiation stage Ogr when 
ra^jrax = 1 in Figure 0(a). In the hard binary stage, if the loss cone in spherical galaxies is depleted, most of the stars that diffuse 
into the loss cone have the energy £^peak = '/'(''peak) Ei^nd the apocenters of the radial orbits, rpeaj;, are labelled as solid triangles. For 
hard BBHs, we have rp^ak 2* rna^(''colli ^h) ^ /ai (/a — 1), which helps to justify the generalization of solution of the Fokker-Planck 
equation in the tidal disruption context to the BBIf systems. 
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5.3.2 Evolution of BBHs with unequal BH masses 

For the BBH with unequal masses, as discussed in § |3.2[ low-angular momentum stars may have been depleted before the 
BBH first becomes hard because the BBH becomes hard at smaller semimajor axes and the total stellar mass in the full loss 
cone decreases with the decrease of the semimajor axes. Figure ^(b)-(d) shows the evolution of the BBH with unequal masses 
(m2/mi < 1). 

As seen from Figure ^b)-(d), roughly, the overall evolution trend of the BBHs with unequal masses looks like those with 
equal masses; but the bottleneck stages of the BBHs with unequal masses can be either in the non-hard binary stage or in 
the hard binary stage, and tend to occur at smaller semimajor axes. The bottleneck timescales of most of the BBHs are still 
longer than the Hubble time. Generally, the bottleneck timescales decrease with decreasing mass ratios in most of the core 
galaxies, but remain unchanged in most of the power-law galaxies. 

Before the BBHs become hard, at large radii (r ~ 10^-10'* pc), the BBH evolution timescales increase with decreasing 
"^12 /mi , but most of them are still smaller than the Hubble time tnubbic = 10" yr (Fig. Ija-d). But when m2/mi <, 0.01, the 
BBH evolution timescales in some of the galaxies exceed the Hubble time tnubbic = 10" yr at r ~ lOkpc. So, in BBHs with 
m-2/mi ^ 0.01, the orbit of the smaller BH would not decay to the galactic center from the outer edge of the galaxy within a 
Hubble time. 

For the BBH with unequal masses, if the removed mass from the galactic core during the non-hard binary stage can be 
ignored and the loss cone is approximately full when the BBH becomes hard, the evolution of the BBH after a = Oh would 
be like that of the BBH with equal masses, which is first controlled by the depletion of the loss cone (when the evolution 
timescale increases sharply) and then controlled by the stellar difi^usion (see flat parts of the BBH evolution timescale curves 
in Fig. ^-d) and gravitational radiation. The bottleneck stage of the BBH with unequal masses shifts to smaller semimajor 
axes because it becomes hard at smaller semimajor axis. Note that the flat part (controlled by stellar diffusion) of the BBH 
evolution curve of each galaxy, if any, has almost the same height in all the panels of Figure |^, which refiects the fact that 
the BH mass ratio does not affect the hardening time (see equation 

For the BBH with unequal masses, if the loss cone is depleted when the BBH first becomes hard [??(£') ~ in equation 
p^ , the following BBH evolution will be directly controlled by the stellar diffusion to the loss cone (see the flat part of the 
BBH evolution curves, mostly for power-law galaxies) and then gravitational radiation, or directly controlled by gravitational 
radiation (no flat part in the BBH evolution curves, mostly for core galaxies). Before the BBH becomes hard, we obtain the 
place adp where the low-angular momentum stars are depleted by setting Mic(adp) — Mmi — Jrm(adp) (see equation ^jj) . 
And we use equation (^6|) (log-linear approximation) to obtain the BBH hardening timescales from Odp to ah because the 
evolution at Oh < a < ffldp is an intermediate process between the evolution controlled by dynamical friction and the evolution 
controlled by three-body interactions with the stars that diffuse from high-angular momentum orbits; then, we use equation 
( prj) to obtain the BBH evolution timescales. As seen from Figure |^(b)-(d), the evolution timescale curves of these BBHs 
increase steeply from a — adp(> flh) and smoothly connect to the timescales in the hard binary stage. If the BBH evolution 
is controlled by stellar diffusion to the loss cone when it just becomes hard, its lifetime is not affected by mass ratios. For 
some of the BBHs, gravitational radiation may become the dominant dissipative force when a > Oh, and in the non-hard 
binary stage their evolution timescales show a peak or maximum (i.e., the bottlenecks) at a = Ogr where tgr(agr) = th(aigr). 
With decreasing BH mass ratio, the merger timescale becomes shorter because the gravitational timescale tgi- (equation [ic| ) 
decreases as a^, though at a given a, tgi is inversely proportional to mass ratio m2/mi (when m2/mi ^ 1). Here, as the BH 
mass ratio decreases, the bottleneck stage shifts to smaller semimajor axes because either ah or adp decreases. 



5.3.3 Sharp change of the BBH evolution timescale at a = Oh 

As seen from Figure |^(a), the BBH evolution timescale increases sharply at a = ah. This sharp increase arises because the 
loss cone is depleted almost as soon as the BBH becomes hard. A similar sharp increase starting at a = Oh also exists in 
Figure |^(b)-(d) for those galaxies whose loss cones are approximately full when the BBH flrst becomes hard. In this subsection, 
we will see that the coincidence between the increase in the evolution timescale and Oh is caused by some simplification in our 
analysis; but even if more realistic treatment is considered, the sharp increase would also start at some place close to a = Oh, 
though not exactly at a = ah, and our conclusions would not be significantly affected. 

Note that in the analysis in § ^ and § ^, before the BBH becomes hard, the velocity dispersion of the low-angular 
momentum stars not removed from the galactic core is assumed to be unaffected by interactions with the BBH; and after the 
BBH becomes hard, each star having had close encounters with the BBH is removed from the galactic core with an energy 
gain. In fact, the heating of the stars is a continuous physical process, i.e., in the non-hard binary stage, the stars should also 
be heated as the BBH hardens even though they are not removed from the galactic core. 

From scattering experiments with the restricted three-body approximation, the BBH energy change due to stellar en- 
counters in the non-hard stage is generally very small [e.g. \AE{a — 10ah)| — 0.03|A_E(a = ah)| — 0.1m«(j^, cf. fig. 1 in 
Quinlan (1996)] and decreases steeply with increasing a (which is consistent with the steep decrease of the stellar removing 
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rate from the core with increasing a at a > Oh in Fig. |l|); only when a is close to ah does the relative change in the BBH energy, 
AE/E, increase steeply towards a constant (independent of E) as the BBH hardens [cf. fig. 1 in Quinlan (1996)]. Hence, in 
the non-hard binary stage, the energy increase of a star occurs through a slow diffusion process. If all the energy loss of the 
BBH is used to heat the low-angular momentum stars, the heating would become significant after Gm\m2/2a ~ 3A/ic(a)(7c/2, 
i.e., after 



5ah 



2mi 



^1/2 


' 0.03M. " 




_Mic(ah)_ 



1/2 



^ 5ah, 



(51) 



where Mic(a) ~ Mic(ah)a/ah for a <C an (cf. equation |l^) and Mic(ah)/M, ~ 0.03-1 for those galaxies whose loss cone are 
approximately full at a = ah in Figure |l|(a). Thus, the heating of the stars becomes significant only at some place close to ah 
and the sharp increase of the BBH evolution timescales would also start not far away from a = ah. The semi-major axis ah 
is also the point where the typical energy gain of a scattered star is just large enough for it to escape from the galactic core. 
Thus, we expect a sharp increase in the evolution timescale near a ~ ah. 



5.3.4 Testing assumptions 

We can use the results obtained above to test some of the assumptions that we have made in our analysis. 

(i) In the hard binary stage, if the loss cone is depleted, the BBH decay is controlled by the stellar diffusion into the loss 
cone. We define Speak as the energy at the peak of stellar diffusion rates F^'^{£,a) (see equation and fpoak is generally near 
or smaller than ■(/'(an), not at fic (a) because it is the factor -Fmax(5) that dominates the shape of F^'^{£,a) in equation ( |2o| ) 
(see also MT). As seen from Figure BI most of the stars contributing to the BBH orbital decay have energy f peak ^ GM,/ fad 
or rpcak ^ max(rcoii, flh) > /ad, where rpoak is defined by fpoak = V'('"poak) and rcoii (see equation |25|) is the collision radius, 
which helps to justify the generalization of solution of the Fokker-Planck equation in the tidal disruption context to the BBH 
systems. Figure ^ also shows that the stellar mass within rpcak is higher than O.lAf. for all the galaxies, which confirms that 
resonant relaxation is not important here. 

(ii) In the hard binary stage, when stellar diffusion dominates the BBH orbital decay, the difl:usion timescale into the loss 
cone at £ ~ fpcak is given by 

Mt(rpcak) Jic(£'poak, /aa) 

tdiff 



/i^(£:,a)d£: J|(£pcak) 

(0.06^ 0.6) (i^) (H^) (^) 



(52) 



which is less than the evolution timescale |a/a|. Thus our use of time-independent solutions to the Fokker-Planck equation is 
justified. 

(iii) With equation (^), we may also obtain the total mass of the stars strongly interacting with the BBHs during the 
hard binary stage, which is less than the core mass as seen from Figure ^ (m2/m\ — 1) and is much less when m2/m\ < 1. 
The removed mass from the galactic core during the non-hard binary stage is also small (~ O.IM,), as seen from Figure ^ 
So, the approximation that we ignore the response of the galactic structure to BBH evolution is reasonable. 



(iv) The center of mass of BBHs is not expected to be located exactly at the galactic center (Bahcall fc Wolf 1976). The 




BBH center of mass should be in equipartition with the stars, which implies that its rms radial excursion in core galaxies is 

(53) 

where o, is the velocity dispersion of the BBH center of mass. The "break radius" rb is roughly equivalent to the core radius 
Tc in core galaxies. Power-law galaxies have no well-defined core, but replacing Vc with the break radius rb in equation dss] ) 
should give an upper limit to the likely motion of the BH. Hence, for both the core galaxies and the power-law galaxies, we 
set rc = rb to obtain the wandering amplitude of the center of mass of BBHs. Wandering of the BBH center of mass might 
enlarge the loss cone, increase the clearing rates and decrease the BBH hardening time (e.g. Quinlan & Hernquist 1997). We 
believe that this effect is generally not important in these spherical galaxies for the following reasons: (1) For all the core 
galaxies, gravitational radiation becomes dominant before their BBH (0.01 ^ m2/m\ ^ 1) semimajor axes decrease to the 
value rran (e.g. Fig. |^). (2) For most of the power-law galaxies, a = rran are located at the flat parts of the evolution curves 
(see Fig. ^) where stellar diffusion dominates the BBH orbital decay, or a — - ?^ran are less than or around the semimajor axes 
where the loss cones are depleted (e.g. for the BBHs with lower mass ratios 0.01 ^ m2/mi < 1), and so the BBH lifetime 
is not sensitive to the wandering of the BBH. In Figure ^a), for only two power-law galaxies, a — rran happens before their 
loss cones are depleted, but their BBH lifetimes are shorter than the Hubble time even without considering the wandering 
of the BBHs. (3) With decreasing m2/m\ ratios (say, 1x12/ m\ 0.001), though a = rran may happen before the loss cone 
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Figure 6. The core mass versus the total mass of stars interacting strongly with BBHs (rjte/mi = 1, in spherical galaxies) during the 
hard binary stage agr < a < Qh (^^^ equation |l2| for T^fintoract^ pjg ^ £qj. Q-jj/agr or Fig. ^ for the BBH evolution). The solid circles 
represent core galaxies and the open circles represent power-law galaxies. For almost all of the galaxies, M«{rij) > jv/i"*'^'''"^' , which helps 
to justify our neglect of the reaction of the core to the inspiral of the BHs. 



is depleted, the dynamical friction timescales at large radii (<; lOkpc) are longer than the Hubble time and the smaller BH 
cannot sink into the center of the galaxy with mi. 



5.4 Axisymmetric galaxies 

Assuming that the galaxies in Table 1 are axisymmetric, we use equations (^-(|ll|) and ( |2^)-(|35| ) to obtain the BBH evolution 
timescales as a function of BBH semimajor axes after the BBH becomes hard (equations [L9|-|2^inay also be used if the galaxies 
are nearly spherical). Before the hard binary stage, the BBH evolution timescales are obtained from the equations in § ^ as in 
spherical galaxies. We will first consider the evolution of BBHs with equal masses (rn^lrrix — 1) in the galaxies with different 
degrees of flattening e (defined in § 3.2.2) (Fig. |^; and then consider the evolution of BBHs with different m2/mi ratios for 
a given value of e(= 0.3) (we estimated in § that the average e in our sample was ~ 0.36) (Fig. 

As seen from Figure |l| or § |5!|, if the BBH has equal BH masses, the removed stellar mass from the galactic core during 
the non-hard binary stage can be ignored and the loss wedge is approximately full when the BBH becomes hard [ri{£) ~ 1 
in equation As seen from Figure ^ (m2/mi = 1), the BBH evolution timescales follow similar trends (small evolution 
timescales at large and small radii, and bottlenecks at intermediate radii) to those in spherical galaxies. The bottleneck occurs 
during the hard binary stage. After the BBHs become hard, the effects of flattening on the hardening timescales depend on 
the value of e. Figure |^ shows that the bottleneck timescales in the most nearly spherical galaxies (e.g. e=0.01, 0.05 and 0.1 
in panels a, b and c) are still longer than the Hubble time, but the bottleneck timescales in most highly flattened galaxies 
(e.g. 21 out of our sample of 30 galaxies in panel d with e — 0.5) become shorter than the Hubble time. If e is small (e.g. 
^ 0.01), for most stars on centrophilic orbits, the characteristic angular momenta Js (see § 3.2.2) is not significantly larger 
than Jic(a) and the evolution reduces to the spherical case (see most of the power-law galaxies or some core galaxies in 
Fig. ^ and the galaxies in Fig. ^). As e increases, the effects of flattening become significant with the increase of the total 
stellar mass in the loss wedge (cf. Fig. As seen from Figure ^b)-(d), the semimajor axes of the BBHs decrease more in 
axisymmetric galaxies than in spherical galaxies before they reach the bottleneck stages, and the bottleneck stages shift to 
smaller semimajor axes as e increases. For many BBHs in nearly spherical or low-e axisymmetric galaxies, (e.g. most of the 
BBHs in power-law galaxies in Fig. |^-c or some BBHs in core galaxies in Fig. ^ whose evolution curves have flat parts at 
intermediate radii ~ 10~'^-lpc), gravitational radiation becomes dominant still after the loss wedge/cone is depleted; and the 
lifetime is controlled by the stellar refilling rate to the loss wedge/cone. As in spherical galaxies, after the loss wedge/cone 
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Figure 7. BBH evolution timescales for axisymmetric galaxies with 1712/ mi = 1 and varying degrees of flattening e. The curves have 
the same meanings as those in Figure ^ In panels (a) and (d), the solid circles (for core galaxies) and the open circles (for power-law 
galaxies) represent the maximum Brownian motion magnitudes of the BBH centers of mass. The evolution curves before the hard binary 
stage are the same as those in Figurej2l The bottleneck is in the hard binary stage. When e is small (e.g. panel a), the BBH evolution 
reduces to the spherical case (cf. Fig. Hi). As e increases, the bottlenecks shift to smaller semimajor axes and many BBH lifetimes become 
shorter. The Brownian motion of BBHs does not affect the BBH lifetime in nearly spherical galaxies (e.g. panel a) or in core ga laxies 
(e.g. panels a and d), but may decrease the BBH lifetime in highly flattened power-law galaxies (e.g. panel d). See details in § t>A. 
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is depleted, the inner parts of most of the axisymmetric galaxies are also in the diffusion regime before the gravitational 
radiation stage begins. In the diffusion regime, the stellar diffusion rate to the loss wedge is also insensitive to a, and hence 
the BBH evolution timescales controlled by the stellar diffusion are nearly constant. As seen from Figure |^a)-(c), the flat 
parts of the evolution curves in Fig. ^ have almost the same heights as those in Figure ^ and so the lifetimes of those BBHs 
are affected little by the flattening effect, which is consistent with the results in MT that the stellar tidal disruption rates by 
central BHs are not significantly changed by the flattening effect. With increasing e, more and more BBH evolution curves 
have no flat parts, i.e., gravitational radiation becomes dominant before the loss wedge is depleted (e.g. BBHs in power-law 
galaxies in Fig. ^ and BBHs in core galaxies in Fig. |^-d); and their merger timescales decrease, because the total stellar 
mass in the loss wedge increases and gravitational radiation becomes dominant at smaller semimajor axes. 

Figure |^ shows the evolution for the BBHs with unequal BH masses when the flattening parameter of the galaxies e is 
0.3. In this case the removed stellar mass from the galactic core during the non-hard binary stage can be ignored and the 
loss wedge is approximately full when the BBHs (with m2/mi ^ O.OI) become hard (see Fig. ^). As seen from Figure ^ 
the bottleneck stages (also in the hard binary stage) shift to smaller semimajor axes as the mass ratio decreases, but not so 
obviously as those in spherical cases (Fig. |^ since the stellar mass in the loss wedges of the galaxies are generally so much 
not to be depleted even when the BBHs reach the gravitational radiation stage. As in spherical galaxies (Fi g. the merger 
timescales decrease with decreasing mass ratios because the gravitational radiation timescale tgi (equation hd) decreases as 
a*. When m2/m\ decreases to O.OI, most of the BBH merger timescales are shorter than the Hubble time (Fig. pn). 



The justifications of our assumptions in § 5.3.4 can be generalized to axisymmetric galaxies (and triaxial cases below). 
After the loss wedge is depleted, most of the stars contributing to the BBH decay have energy fpcak ^ GM,/faa or rpcak ^ 
max(rcoii, a-h) > fai; the stellar diffusion timescales to the loss wedges are smaller than the BBH hardening timescales; and the 
removed stellar mass from the core is also less than the total core mass. If e is small (e.g. e = 0.01 in Fig. ^), the wandering 
of BBHs is not important, just as in spherical galaxies. When e is large (e.g. e = 0.3 or 0.5 in Fig. ^ or ^d), the wandering of 
BBHs is also not important in core galaxies. In a few of the power-law galaxies whose lifetimes are longer than the Hubble 
time, a = rran happens before the loss wedges are depleted and their BBH evolution timescales at a = rran are shorter than 
the Hubble time, so the wandering of BBHs might decrease the BBH lifetime in those galaxies. 



5.5 Triaxial galaxies 



Using the equations in § 3.2 and § ^, we obtain the BBH evolution timescales if the galaxies in Table I are triaxial. Here, we 
mainly consider the BBH with equal masses (m2/m\ — 1) to see how the evolution depends on the triaxiality e. The trends 



of the BBH evolution with different 1112/ m\ are similar to those in axisymmetric galaxies (see Fig. g or § 5.4). Noti ng tha t 
triaxial galaxies containing central BHs are likely to evolve secularly toward axisymmetry within a time X'''^^"^ (see § 3.2.3 ), 
we will first ignore the evolution toward axisymmetry by setting X'''^^"^ = 00 and see the maximum effect of the triaxiality 
on the BBH evolution. The results are shown in Figure H. Then, we will set finite y*"^^"" to see the effect of triaxiality on the 
BBH evolution. The age of galactic triaxiality T^gJ, is assumed to be zero at a = Oh. 

Figure H shows the BBH evolution in triaxial galaxies without including secular evolution towards axisymmetry. As seen 
from Figure H, the BBH evolution curves in triaxial galaxies follow similar trends to those in spherical or axisymmetric galaxies 
and the effects of triaxiality on the evolution timescales also depend on the value of e as in axisymmetric galaxies. If e is 
small, the characteristic angular momentum Js is not significantly larger than Jic and the BBH evolution would be similar 
to that in spherical or axisymmetric cases. For example, when e = O.OI (Fig. ^a), most of the BBH evolution timescales of 
the power-law galaxies are like those in the spherical or axisymmetric galaxies and the bottleneck stages are controlled by 
the stellar diffusion to the loss cones/wedges. When e is large (<; 0.05, see Fig. ^s-d), the difference between the loss region 
and the loss cone/wedge becomes significant. The loss regions in all the galaxies cannot be depleted before the gravitational 
radiation stage begins because there are many more stars in the loss regions J < Js than in the loss cones/wedges. Most of 
the BBH lifetimes in Figure |^(b) (e = 0.05) are comparable to or shorter than the Hubble time and almost all the BBHs in 
Figure |(c)-(d) (e = 0.1 and 0.5) can merge within a Hubble time. As in axisymmetric galaxies, only the lifetime of BBHs in 
power-law triaxial galaxies with large e (e <; 0.05 in panels b-d) might be decreased by their wandering motion. 

After including secular evolution towards axisymmetry, the effect of triaxiality on the BBH evolution depends on the 
value of T*"'"'"". For example, if r*'''"'' = W^P{£), after the BBHs become hard, they first follow similar evolution to those 
in Figure then with evolution towards axisymmetry, they evolve as those in axisymmetric galaxies (Fig. ^ ; and the BBHs 
have similar lifetimes to those in axisymmetric galaxies. If is smaller, the BBH evolution curves would be close to 

those in axisymmetric galaxies at earlier time; and if is larger, they would be more like those in Figure ^ Actually, 

if y*'^''"'' = 10^P(f ), the BBH evolution curves in triaxial galaxies with large e (e.g. ^ 0.05) are quite similar to those in 
Figure |(b)-(d). 

In the discussion below on BBHs in triaxial galaxies, we will always use the results in Figure ^, which gives the maximum 
effect of triajciahty on BBH evolution. 
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Figure 8. BBH evolution timescales for axisymmetric galaxies with e = 0.3 and varying mass ratios. The curves have the same meanings 
as those in Figure |2| In panel (a), the solid circles (for core galaxies) and the open circles (for power-law galaxies) represent the maximum 
Brownian motion magnitudes of the BBH centers of mass. As the mass ratio decreases, the bottlenecks shift to slightly smaller radii and 
the merger timescales become shorter. Brownian motion of BBHs does not affect the BBH lifetime in core galaxies, but may decrease 
the BBH lifetime in power-law galaxies. 
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Figure 9. BBH evolution timescales for triaxial galaxies with m2/m\ = 1 and varying degrees of triaxiality e. Secular evolution towards 
axisymmetry is not included. The curves have the same meanings as those in Figure bl The solid circles (for core galaxies) and the 
open circles (for power-law galaxies) represent the maximum Brownian motion magnitudes of the BBH centers of mass. Comparing with 
spherical galaxies and axisymmetric galaxies (Fig. |^ and |7|, the draining of the loss regions in triaxial galaxies can significantly decrease 
the BBH lifetime. Almost all the BBHs can merge within a Hubble time when e is large (panels c and d). The lifetime of BBHs in the 
power-law triaxial galaxies with large e (e.g. ^ 0.05 in panels b— d) may be decreased by their Brownian motion. 
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5.6 Generalization to a distribution of stellar masses and radii 

In the calculations above, we assume a single stellar mass of IM© and a single stellar radius of IKq. In this subsection, we 
will see that generalizing to a distribution of stellar masses and radii does not significantly affect our results. 

First, generalizing to a range of stellar mass does not affect the dynamical friction timescales (see equations ^^-^^ , which 
is affected by the mass density, not the mass of a single star. The value of A may be sensitive to the stellar mass, but the 
friction timescale depends only logarithmically on A. 

Second, in the hard binary stage, the BBH hardening timescale is controlled by the mass clearing rate from the loss region 
(see equation |^ . If the loss region is not depleted, generalizing to a range of stellar mass does not affect the BBH hardening 
timescale since the total mass cleared from the loss region per unit time is not changed. Only when the evolution is controlled 
by stellar diffusion to the loss region is the BBH evolution timescale affected through the variation of two-body relaxation 
diffusion coefflcient with stellar mass (see equation ^) . According to equation ( A6) in MT, the diffusion coefficient of two- 
body relaxation at constant mass density fl cx m^n(m*)dm,, where n(m,) is proportional to the probability of finding a 
star with mass m* m* -f dm, and n(m,)dm, = M0. Thus, if we use a Salpeter mass function n(mt) oc mj^'^^ with 
limits mi — 0.08 M© and m2 = 1 M©, the diffusion coefficient is reduced by a factor of 0.31 compared to a population of stars 
of mass 1 M0. Therefore, when the stellar diffusion to the loss cone or loss wedge controls the BBH evolution timescales, the 
evolution timescales will be only increased by a factor of three or so, which does not affect our qualitative conclusions. 

Third, if the loss region is depleted in the non-hard binary stage, the evolution of the BBH after the depletion and before 
the hard binary stage is controlled by both dynamical friction from distant stars and three-body interactions with low-angular 
momentum stars diffused from high-angular momentum orbits, which further reduces the effects of the dependence of the 
diffusion coefficient on the stellar mass distribution. 

Finally, the range of stellar radii has relatively little effect on the collision radius radius rcoii (equation ^). Thus, we still 
have Tcoii <S rpoak in Figure ^ and the BBH evolution timescale is not affected much (see § 5.2). 



Since generalizing to a distribution of stellar masses and radii does not significantly affect the BBH evolution timescales. 



we will continue to use the results with a single stellar mass and radius (obtained in § 5.3 -§ 5.5) to study the properties of 
surviving BBHs in § M. 



6 PROPERTIES OF SURVIVING BBHS 



As seen from § ^, BBHs can survive for a Hubble time in spherical galaxies, axisymmetric galaxies or some galaxies with small 
triaxiality (e ^ 0.05). Whether or not some of the various astronomical phenomena that have been associated with BBHs [e.g. 
double nuc l ei; double-peaked Balmer lines ([Caskell 1996 ); quasi-periodic radio, optical. X-ray or 7-ray variability ( Sillanpaa 



et al. 198S 



Rieger fc Mannheim 2000) etc.] can be observed depends on the current orbital properties of the BBHs, such as 



their semimajor axis a, circular speed Vc — ^ GM,/a, orbital period Porb = 2-Ka/vc, etc. 

Figure ^ shows the relation between the galactic velocity dispersions ctc and the BBH orbital parameters at the point 
where \a/a\ = 10^" yr as a decreases for the first time in Figure ^ (spherical galaxies). As seen from Figure [ic|, the BBH orbital 
parameters depend on galactic velocity dispersions and BH mass ratios. The BBHs in most of the galaxies with CTc ^ 90 kms~^ 
have merged. The semimajor axes and orbital periods of surviving BBHs are generally in the range 10~^-10pc and 10-10^ yr; 
and they are generally larger in high-dispersion galaxies (core galaxies with high central BH masses) than in low-dispersion 
galaxies (power-law galaxies with low central BH masses), and larger for BBHs with equal BH masses than for BBHs with 
unequal masses. The orbital velocities of surviving BBHs are generally in the range 10^-10* kms~^; and they are generally 
larger in high-dispersion galaxies than in low-dispersion galaxies and larger for BBHs with unequal masses than for BBHs 
with equal masses. The surviving BBH properties estimated above may help identify appropriate methods to probe BBHs in 
different spherical galaxies: 

(i) Double nuclei associated with BBHs should be easier to observe in high-dispersion or luminous galaxies than in low- 
dispersion or faint galaxies because BBH semimajor axes are much larger in high-dispersion galaxies. The BBH sample 
identified from double nuclei would have a bias toward high-m2/mi (say, <; 0.1) BBHs for at least three reasons: (1) the two 
components have comparable brightness if the stellar or non-stellar luminosity is correlated with the BH mass, (2) the BBH 
semimajor axes are generally larger if the BH masses are equal. (3) BBHs with very low-m2/mi ratios (say, 0.001) may not 
exist in the galactic center since the smaller BH may not sink to the galactic center in a Hubble time (see Fig. ^) (here we 
only consider the BBH formed by mergers of galaxies, and we do not consider the possible low-mass BH formation caused by 
some mechanism at galactic centers). 

(ii) Broad line regions (BLRs) associated with the two components of a BBH in AGNs may produce double-peaked emission 
lines (e.g. Balmer lines) varying with timescales of the BBH orbital period, though the double peaks are not easily recognized 
for Vc ^ (Jblr where (tblr is the velocity dispersion of the clouds in BLRs (typically 3000-5000 kms~^). Though the BBH 
orbital period can be as short as 10 yr for mijrax = 0.01 (Fig. |lOp), for such an extreme mass ratio one component of double 
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Figure 10. Estimated observational properties of surviving BBHs in spherical galaxies with varying mass ratios. The panels show the 
relations between the galactic velocity dispersion Cc and the BBH semimajor axes a, circular speeds Vc = ^ GM, / a and orbital periods 
^orb = 2na/vc at the point where \a/a\ = 10^" years for the first time in Figure |2| The solid circles represent the core galaxies and the 
open circles represent the power-law galaxies. The solid triangles (core galaxies) and open triangles (power-law galaxies) located at the 
bottom of each panel represent the dispersion of the galaxies in which the BBH evolution timescales are less than the Hubble time 10^" 
years. BBHs in most of the low-dispersion (ctc ^ 90kms~^) galaxies have merged within a Hubble time and BBHs are more likely to 
survive in the galaxies with ^ 90kms~^ (similarly in Fig. llHl3 below). The semimajor axes and orbital periods of surviving BBHs 
are larger in high-dispersion galaxies than in low-dispersion galaxies (similarly in Fig. below) and larger for BBHs with equal 

masses than for BBHs with unequal masses (similarly in Fig. 12 below). The circular speeds are larger in high-dispersion galaxies than 
in low-dispersion galaxies and larger for BBHs with unequal masses than for BBHs with equal masses (similarly in Fig. ^ below). Note 
that BBHs with very low mass ratios (say, ^ 0.001) arising from mergers of ga laxies would not exist at galactic centers because the 
smaller BH cannot sink into the galactic center within a Hubble time (cf. § 5.3.1). 
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peaks would be too weak to be discernible compared to the other strong component. The double peaks from the BLRs of 
high-m2/mi BBRs would have comparable strengths, but it is hard to detect the line variability within a reasonable time for 
high-m2/mi BBHs since Porb is at least 10^-10'' yr (see the case m2/m\ = 0.1 in Fig. [lo|b) even in low-dispersion galaxies 
(o"c ^ 150 kms"^). So, double-peaked emission lines from BLRs are unlikely to show periodic behavior due to the long orbital 
period in BBH galaxies. For example, the quasar 3C390.3 was introduced as a BBH candidate based on data from only ~ 20 
years, which is much less than the estimated BBH orbital period Porb ~ 300 yr. The radial velocity curve of ~ 20 yr for one 
peak of emission line H/3 in 3C390.3 is consistent with the expectation from the orbital motion of a BBH ( [Gaskell 1996|); but 



the results from continuous observations of ^ 10 more years deviate from the expectation from a BBH model (Eracleous et 
al. 19f|). 

(iii) Some periodic behavior observed in the radio, optical. X-ray or 7-ray lightcurves is possibly related to a BBH creating 
jet(s) aligned nearly along the line of sight or interacting with disk(s). The observed period is usually identified with the 
orbital period of the BBH or a fraction of the orbital period. For example, if a relativistic jet emerges from the less massive 
BH and is aligned nearly along the line of sight, the flux of the X/7-ray radiation emitted from the jet may vary with the 
periodic change of the Doppler-boosting factor (due to the slight change of the jet inclination angle with the BH orbital 



motion) (Rieger & Mannheim 2000). Because of the same mechanism in the phenomenon of "superluminal motion" (i.e. the 
observed timescale is shorter than the intrinsic timescale, see Rees 1966), the observed X/7-ray variability timescale is shorter 
than the orbital period [say, ~ 10~^Porb in Rieger & Mannheim (2000)] which can be identified within a much shorter time 
(~ 0.1-1 yr) in low-dispersion galaxies with m2/m\ ~ 0.01 (Porb ~ 10-lOOyr when Oc ^ 150kms~^ in Fig. p^). It is hard 
to use other variability phenomena with timescales ~ Porb to search for BBHs within a short time (say, < 50 yr) except for 
low-m2/mi (~ 0.01) ones in low-dispersion galaxies (ctc ^ 150kms~^). 



In axisymmetric or weakly triaxial (e ^ 0.05) galaxies (Fig. [11|-|13|), the relations between the orbital parameters of 
surviving BBHs and the properties of their host galaxy follow similar trends to the relations in spherical galaxies (Fig. |lo[ ), 
though the normalization may be different. The orbital properties of surviving BBHs are in the similar range as those in 
spherical galaxies, and they depend on the fiattening or triaxiality parameter e, as well as the galactic velocity dispersion 
and BH mass ratio. As seen from Figure (BBHs with mass ratio m2/m\ = 1 in axisymmetric galaxies), with increasing 
flattening parameter e, BBH semimajor axes decrease (the upper limit is ~ 10 pc for e = 0.01 and ~ Ipc for e = 0.5), which 
may require higher resolution to resolve double nuclei; BBH orbital periods decrease (~ 10^-10^ yr for e — 0.01 and ~ 50- 
10^ yr for e — 0.5), which may help in the identiflcation of periodicities in X/7-ray variability within a short time (~ 0.5-10 yr 
if the timescale of periodicities ~ 10~^Porb), but the orbital periods are still so long that the periodic variability of double- 
peaked emission lines from BLRs will be difficult to identify. As seen from Figure |l^ (flattened galaxies with t = 0.3), with 
decreasing m^lmx ratios, the BBH semimajor axes and the orbital periods decrease as in spherical galaxies. When 1712/ mi 
decreases to 0.01 (Fig. [l^), many BBHs (22 out of our sample of 30 galaxies) have merged within a Hubble time; and most 
of the surviving BBHs in low-dispersion galaxies {a^ ^ 150kms~^) have an orbital period ~ 10-30 yr and could possibly be 
found by identifying some variability phenomenon associated with the orbital motion. It is difficult to observe BBHs even 
with equal masses in more triaxial galaxies (e <; 0.1) because most of them have merged within a Hubble time (Fig. |l^ and 
d). In addition, the random motion of the center of mass of BBHs might affect their observational properties in low-dispersion 
axisymmetric galaxies (with e ^ 0.1) or low-dispersion triaxial galaxies (with e ~ 0.05). On one hand, the random motion may 
decrease their semimajor axes and orbital periods, which would make variability phenomena a more efficient tool to search 
for BBH candidates; on the other hand, it may also make more BBHs merge within a Hubble time and decrease the number 
of surviving BBHs, which would decrease the probability of identifying BBHs by variability phenomena. 

Note that the double-peaked emission lines from BLRs and radio, optical, X/7-ray variability that we have discussed 
above are phenomena seen in the relatively rare and distant AGNs, not the common inactive nearby galaxies where BHs are 
known to exist. 

If future observations do not reveal any phenomena associated with BBHs in AGNs, that might suggest that: (i) the 
inner stellar distribution in AGNs is different from normal nearby early-type galaxies, (ii) gas in AGNs plays an important 
role in BBH orbital decay (e.g. Gould & Rix 2000), or (iii) nuclear activity switches on soon after a galaxy merger and the 
lifetime of nuclear activity is much less than the Hubble time, so the surviving BBHs in AGNs would not have had time to 
spiral together before the nuclear activity ceases. 

To see if double nuclei associated with BBHs formed by mergers of galaxies can be observed with current telescope 
resolution, we show the predicted BBH semimajor axes m arcsec versus the galactic velocity dispersion CTo for the nearby 
inactive galaxies (Table 1) in Figure |l^. The BBHs are assumed to have equal BH masses, which have larger semimajor 
axes than those with unequal masses. As seen from Figure spherical or nearly spherical and high-dispersion galaxies have 
surviving BBHs with the largest semimajor axes, which are just within the HST resolution, 0.1". 
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Figure 11. Estimated observational properties of surviving BBHs in axisymmetric galaxies with ni2/nii = 1 and varying degrees of 
flattening e (see Fig. ^). The symbols have the same meanings as those in Fig. The relations between the surviving BBH orbital 
properties and the galactic dispersion follow similar trends to the relations in spherical galaxies (Fig. . Most BBHs in highly flattened 
galaxies can merge within a Hubble time (panel d). The semimajor axes and the orbital periods of surviving BBHs are larger in nearly 
spherical galaxies than in highly flattened galaxies; while the circular speeds are larger in highly flattened galaxies than in nearly spherical 
galaxies (similarly in Fig. [13 below). In power-law or low-dispersion galaxies with large e (e.g. panel c or d), Brownian motion of BBHs 
may further decrease the semimajor axis (at the same time, decrease the orbital period and increase the circular speed) shown in the 
flgure or make the BBHs merge within a Hubble time (similarly in Fig. 12H14 below). 
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7 DISCUSSION AND CONCLUSIONS 

We have studied the evolution of BBHs in realistic galaxy models obtained from a sample of nearby early-type galaxies. 
We calculated the BBH evolution timescales due to stellar interactions, taking into account the refilling of the loss cone by 
two-body relaxation and by tidal forces in non-spherical galaxies. The evolution of BBHs depends on BH mass ratio and host 
galaxy type. BBHs with low mass ratios (say, m2/m\ <, 0.001) are only rarely formed by mergers of galaxies because the 
dynamical friction timescale is too long for the smaller BH to sink into the galactic center within a Hubble time (cf. Fig. |^. 
In spherical galaxies, if the BBH has equal BH masses, low-angular momentum stars are generally depleted quickly after the 
BBH becomes hard and thereafter the BBH lifetime is controlled by the stellar diffusion rates into the loss cone (see Figs. 
^ and^; if the BBH has unequal masses (so that the BBH becomes hard at smaller semimajor axis), low-angular momentum 
stars can be depleted before the BBH becomes hard, and possibly gravitational radiation can dominate the BBH orbital decay 
in the non-hard binary stage and the BBH lifetime decreases (see Figs. ^ and |^-d) . Low-dispersion galaxies have smaller 
central BHs and higher central densities, and hence their BBHs generally have shorter hardening timescales (associated with 
interactions with stars) in the hard binary stage (see equation^ and shorter lifetimes (see Fig. In axisymmetric and triaxial 
galaxies, many stars can precess onto the low-angular momentum orbits by tidal forces and decrease the lifetime of BBHs (see 
Figs, ^-f, l^-j^- The wandering of the BBHs is generally not important, but might decrease the lifetime of BBHs in some 
axisymmetric or triaxial galaxies with low velocity dispersion (see Figs. |^-^. Note that in the study, stars escaping from the 
core are considered as removed from the loss cone, which is usually a plausible assumption; but in some circumstances, even 
a small fraction of stars returning to the loss cone could enhance the BBH decay rate, which deserves further investigation. 

Most of BBHs in the galaxies with velocity dispersion ctc ^ 90 kms^^ have merger timescales shorter than a Hubble time 
(see Fig. |ll[-p^, and so do BBHs in highly flattened (e ^ 0.5) or moderately triaxial (e ^ 0.05) galaxies (see Figs. pT| d and 
p^-d). BBHs with low-m2/mi (say, ~ 0.01) ratios, which have shorter lifetime than the BBHs with equal masses, can also 
merge within a Hubble time in axisymmetric or triaxial galaxies (see Fig. p^). Surviving BBHs are generally in the galaxies 
with velocity dispersion CTc ^ 90kms~^. Spherical galaxies, axisymmetric or weakly triaxial galaxies are all likely to have 
surviving BBHs, especially have those with equal masses (see Figs. |10|-{L2| and |l3|a-b). We also estimated the observational 



properties of surviving BBHs. The semimajor axes and orbital periods of surviving BBHs are generally in the range 10"^-10pc 
and 10-10^ yr; and they are generally larger in high-dispersion galaxies than in low-dispersion galaxies, larger in spherical 
galaxies than in non-spherical galaxies, and larger for BBHs with equal BH masses than for BBHs with unequal masses. The 
orbital velocities of surviving BBHs are generally in the range 10^-10* kms~^; and they are generally larger in high-dispersion 
galaxies than in low-dispersion galaxies, larger in non-spherical galaxies than in spherical/nearly spherical galaxies and larger 
for BBHs with unequal masses than for BBHs with equal masses. 

In short, our study shows that BBHs arising from galaxy mergers are likely to have merged in low-dispersion or highly 
flattened/triaxial galaxies, but should still survive in spherical/nearly spherical and high-dispersion galaxies. The study in 
this paper may help to further explore the merger history of massive BHs. 

The BBH mergers driven by stellar dynamics are also simply discussed by Gould & Rix (2000). They conclude that 
mergers driven by stellar dynamics are nearly impossible in a Hubble time and present a mechanism by which gas can drive 
BH mergers. However, their conclusion is based on much simpler galaxy models and dynamics than those used in this paper. 

We have also discussed methods to detect surviving BBHs. For inactive galaxies, currently, the only method to probe 
BBHs is to search for double nuclei at galactic centers. It would be easier to find double nuclei associated with BBHs in 
luminous (or high-dispersion) nearly spherical galaxies and the BBHs to be found would have a bias toward equal BH masses. 
For active galaxies, we can also use periodic variability phenomena to search for BBHs. Within a short time (e.g. 1 month-1 
yr), it would be easier to find X/7-ray periodic variability associated with the BBH orbital motion in faint or low-dispersion 
(e.g. (Je ^ 150kms~^, but <: 90kms^^ since most BBHs have merged in the galaxies with ^ 90kms^^) galaxies than 
in luminous or high-dispersion (e.g. CTc <: 150kms~^) galaxies. Here, the variability timescale is usually a fraction of the 
BBH orbital period (e.g. 0.01 times of the BBH orbital period), and the BBHs to be found would have a bias toward low 
7712/7711 (e.g. ~ 0.01) ratios. By other variability phenomena with timescales of the BBH orbital period, it is also easier to 
find the low-77i2/7ni BBHs in faint or low- dispersion (ctc ^ 150kms~^, but ^ 90kms~^) galaxies, but it would take at least 
~10-100yr. It is hard to find BBHs by double-peaked emission lines from BLRs either because it takes a long time (^ 50 yr) 
to identify the periodic variability of double peaks with timescales of BBH orbital periods or because the strength of one 
component is too weak to separate it from the other one. 

The upper limit of the semimajor axes of surviving BBHs for all of the galaxies in our study is ~ 10 pc, which is just 
within the HST resolution, 0.1", or several pc for typical galaxies in the sample used in this paper. None of these nearby 
common inactive galaxies has shown signs of surviving BBHs (close double nuclei with separation 0.1") at their centers. It is 
not yet clear whether the absence of detected BBHs in nearby galaxies presents a problem, and if so whether the problem is 
with the calculations in this paper, the hierarchical structure formation model, or the observational detection techniques. 
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Table 1. Galaxy Sample. The second column gives the Hubble types of the galaxies: 1=E, 2=E/S0 or SO, 3=Sa or Sb, 4=dE or 
dSph. The ri, and /^b are the break radius and break surface brightness (corrected for Galactic extinction, in V mag arcsec"'^) 
in the Nuker law (equation [Zt] ); and a, 13 and 7 are the exponents describing the sharpness of the break, the outer and inner 
slopes of the Nuker law. The Vc is the half-light radius and Ty is the mass-to-light ratio in the V band. M, is the central black 
hole mass and CTg is the luminosity- weighted line-of-sight velocity dispersion inside r,, . All the parameters except M, and o-g are 
adapted from Faber et al. (1997) where Hq = 80kms~^ Mpc~^. The as comes from the calculation based on a spherical and 
isotropic model in this paper and Mm is obtained from the correlation of BH mass with galactic velocity dispersion (equation 
|5o| ). The 9 galaxies with "_DF < 0" are deleted from our study because their distribution functions obtained from the Eddington 
formula are negative (equation ^8[ ); and so are the last two galaxies with 7 > 1 because their Ce are significantly affected by 



the gravitational potential of central BHs (see § 5.1) 
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Figure 12. Estimated observational properties of surviving BBHs in axisymmetric galaxies with e = 0.3 and varying mass ratios (see 
Fig. ^). The symbols have the same meanings as those in Fig. The relations between the surviving BBH orbital properties and mass 
ratios follow similar trends to those in spherical galaxies (FigTrd). With decreasing BH mass ratios, more and more BBHs can merge 
within a Hubble time. 
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Figure 13. Estimated observational properties of surviving BBHs in triaxial galaxies with m^lmx = 1 and varying degrees of triaxialitv 
e (see Fig. ^) and without including secular evolution towards axisymmetry. The symbols have the same meanings as those in Figure [lo[ 
Except for some BBHs in weakly triaxial galaxies (e.g. e = 0.01 or 0.05 in panel a or b), most BBHs in triaxial galaxies can merge within 
a Hubble time. 
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Figure 14. Estimated semimajor axes in arcsec of surviving BBHs (with m2/m\ = 1) versus galactic velocity dispersion (Tc. The 
symbols have the same meanings as in Figure |l^. The galaxies are assumed to be spherical in panel (a), axisymmetric in panels (b)— (f) 
and triaxial (without including secular evolution towards axisymmetry) in panels (g)-(j). The upper limit of the BBH semimajor axes is 
just within the i?5T resolution, 0.1". 



© 0000 RAS, MNRAS 000, 000-000 



